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The effects of a displacive structural phase transition on sUding friction are in principle accessible 
to nanoscale tools such as the Atomic Force Microscopy, yet they are still surprisingly unexplored. 
We present model simulations demonstrating and clarifying the mechanism and potential impact of 
these effects. A structural order parameter inside the material will yield a contribution to stick-slip 
friction that is nonmonotonic as temperature crosses the phase transition, peaking at the critical 
Tc where critical fluctuations are strongest, and the sliding-induced order parameter local flips from 
one value to another more numerous. Accordingly, the friction below Tc is larger when the order 
parameter orientation is such that flips are more effectively triggered by the slider. The observability 
of these effects and their use for friction control are discussed, for future application to sliding on 
the surface of and ferro- or antiferro-distortive materials. 



Understanding and controlling nanoscale friction are 
among the top priorities in nanoscience and technology, 
where moving elements are increasingly important. Un- 
earthing mechanisms capable of altering dry friction be- 
tween solids, to be employed in addition to traditional 
means such as lubrication, tuning of load, temperature 
and speed [T], is of great interest, in physics as well as for 
potential applications. The idea is to replace a "dead" 
substrate, whose role is purely passive, with a "live" ma- 
terial hosting a phase transition. Early surface science 
work provided an indirect hunch, in the form of a pre- 
dicted drop of the two-dimensional diffusion coefficient D 
for a brownian adsorbate particle, caused by an underly- 
ing surface phase transition [2]-[4]. In that case, Einstein's 
relation Drf = ksT implies the prediction of a peak in the 
particle's viscous friction r] at the surface critical temper- 
ature T^'^'^^ . Although this linear response mechanism is 
not realistically relevant to dry friction, which is domi- 
nated by stick-slip and strong nonlinearities [1 , that ex- 
ample is nonetheless suggestive of a frictional anomaly 
near a bulk phase transition of the substrate. Experi- 
mentally, a spectacular anomaly, elegant albeit confined 
to a very special case of electronic friction, is the critical 
drop of frictional dissipation (noncontact and contact) 
reported at the superconducting-normal transition of a 
metal substrate [5, 6 . 

In some ferroelectrics such as TGS (triglycine sul- 
fate) [7] and to a lesser extent BaTi03 [8|, AFM topog- 
raphy and friction have shown contrast between surface 
domains. However, the critical temperature dependence 
of sliding friction has not been explored; and these sys- 
tems are pestered with structural and electrostatic com- 
plications which one may wish to avoid at this first, more 
fundamental level. At that level, the basic questions are 
(i) what is the frictional coupling mechanism between tip 
motion and a substrate structural phase transition; (ii) 
what is the distinguishing element of the phase-transition 
related frictional contribution relative to the background 
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FIG. 1: Sketch of the 2D model system, the zoom shows 
the on-site potential symmetries, with the hexagon vertices 
representing the six valleys of the on-site potential, displayed 
below. 



friction; and (iii) could one hope to achieve some friction 
control through external fields that influence the sub- 
strate order parameter? 

In order to explore these questions theoretically ahead 
of future experiments, we resort to a model study. Specif- 
ically, we carry out Molecular Dynamics (MD) classical 
simulations for the stick-slip dry friction of an idealized 
point-like AFM tip over a two-dimensional (2D) model 
crystal substrate whose bulk undergoes a weak, con- 
tinuous structural phase transition. Our chosen model 
will be one of a general "displacive" type, a category in 
principle describing popular systems such as some ferro- 
electrics [Sllig. 

Fig. [T] sketches our model system, a 2D close-packed 
atomic lattice (the substrate) of classical particles (the 
atoms), and a point slider (the tip) pulled over the 
substrate edge (the surface) through a spring (the can- 
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tilever), as in a classic Tomlinson model \V . The sub- 
strate atoms (of mass m) are held together at an av- 
erage distance a by an interatomic pair potential Uij^ 
locally similar to a Lennard Jones (LJ) potential, Uij = 
—U + (^{\ri — Tj \ — a)^ + Pil'^i ~ \ ~ approxi- 
mation which avoids complications including thermal ex- 
pansion. Parameters a and /3 are obtained by fitting 
a LJ potential of depth U and radius a ([/, a and the 
mass m, define our set of "natural" units). In an in- 
terval [—0.1,0.1] centered on the minimum at a, the 
fit yields a = 2S.32U/a'^ and /3 = 7S4..35U/a^. In 
addition to the interatomic interaction Uij, each atom 
is subject to a six-valley on-site potential Ui with the 
same symmetry as the lattice (see inset in Fig. [T]) cho- 
sen such as to favor a small uniform distortion of all 
sites towards one of six equivalent valleys A = 1,...,6 

U. = Um- (3- -4|).f + ^uf. Here 

i^f = + zf is the displacement magnitude of the i-th 
atom from the site center, ao is the distance between the 
minima and the center, and Um and Um are the height 
and the depth of maximum and minima {Um = 0.1/7, 
Um = —0.1/7 and ao = 0.05a). At low temperature T <C 
\Um — /7m 1 7 the ideal 2D lattice is characterized by the 
displacive vector order parameter S{T) = (u) = (r — tq), 
measuring the average distortion u of atoms from the 
central triangular lattice positions tq. At T = 0, the 
substrate minimal energy state is at |^(0)| = ao, with all 
the atoms in the same valley A, and ao/a <C 1. MD simu- 
lations of this model are carried out using 2D rectangular 
cells of large but finite thickness Lz = Nza\/3/2 along z 
and length = N^a (typically A^^; = 40, = 40), first 
in thermal equilibrium, with bulk-like periodic boundary 
conditions (PBC) applied along both x and z; then out 
of equilibrium, with PBC along x alone, the tip sliding 
taking place on one of two edges, in frictional simula- 
tions. We stress here our intent to mimic qualitatively 
the behaviour of a real 3D system with a continuous 
phase transition, with the order parameter coupled to 
the slider's motion [22 . Models similar to the present 
one are routinely and successfully employed in the de- 
scription of the displacive structural phase transitions of 
many systems, notably the very well known ferroelectric 
and distortive ones in the perovskites [9 . In bulk simu- 
lations which we carry out first (details in Supplemental 
Material) the substrate equilibrium structural transition 
is identified by the vanishing (near ksTc = 0.075/7) of the 
order parameter S{T), the peak of susceptibility compo- 
nents Xa/3 = -{{ua,u^) - (ua) (ufs)) / K bT and, slightly 
shifted due to finite size effects, the specific heat peak 
Cv = {{E'^) - {E)'^)/{KbT'^N^Nz), where E is the inter- 
nal energy. Below Tc, symmetry between the six valleys is 
broken and one of the six prevails. Just above Tc, sym- 
metry is thermally restored and the distortion of each 
site, though instantaneously still present, is randomly 
distributed between all six valleys. Near Tc the system 



develops long correlations comparable with the simula- 
tion cell size, and its dynamics becomes correspondingly 
slow, as expected at a second order phase transition. In 
subsequent frictional simulations, the PBC along z are 
removed generating two free surfaces. The point "tip" 
of mass M = 500m and coordinate (X(t),Z(t)), inter- 
acting with the substrate atoms via a LJ potential of 
depth V ~ 0.6/7, is dragged over one surface at con- 
stant velocity vq = b ■ 10~^ ^/U/m through a spring 
of constant k — bU/o?, representing the cantilever lat- 
eral stiffness. A load force Fq (typically of order 1 in 
units of U/a) is applied along z to press the tip onto 
the substrate. The overall equations of motion mX = 
- k{X{t) - vot) and mZ = - Fq are inte- 

grated with a velocity- Verlet algorithm with a time step 
At = 5 • 10~^ Y^ma^//7. The frictional Joule heat is re- 
moved by a standard Langevin thermostat endowed with 
a viscous term —m^r and a corresponding random noise, 
both attached to the slab bottom layer only (7 = 10 
optimally mimics a semi-infinite substrate at an equilib- 
rium bath temperature T [H]). Simulation times are long 
enough for meaningful averages (no less than 50 stick-slip 
events required) , but short enough to avert undesired or- 
der parameter destruction due either to small size, or to 
the nucleation of defects heralding the onset of a BKT 
state. The spring elongation Fx{t) = —k{X{t) — v^t) 
measures the instantaneous friction force. Parameters of 
each frictional simulation are thus a) the temperature T; 
b) the overall substrate order parameter valley A=l,...,6 
for T < Tc] the load Fq] d) the average tip sliding 
speed vq] e) the tip effective parameters such as mass M 
and stiffness k. Among these parameters, velocity is the 
least critical since stick-slip is known to yields a nearly 
speed independent friction coefficient [12 . To cut com- 
putational costs, we generally adopt a rather large speed 
{vq = 0.005^yU/m) - except when good quality stick- 
slip details are needed, requiring slower motion. Two 
tip mass values, M = 500m (results shown here) or 5m 
(shown in Supplemental Material), amounting to a factor 
10 in cantilever frequency ^Jk|M again gave rather sim- 
ilar results. The dependence of friction upon the other 
parameters will be described next. 

Primarily, we examine the lateral spring force Fx{t), 
which exhibits for low speed a classic monoatomic stick- 
slip behaviour (see Supplemental Material), close to that 
observed in realistic AFM nanofriction [12] |^ . Its time 
average measures the dynamic friction and the corre- 
sponding friction coefficient /i = {Fx{t))/Fo. The over- 
all behavior obtained for friction coefficient versus tem- 
perature, order parameter valley, and load - our main 
result - is presented in Fig. |2] First (and not sur- 
prising), the phase-transition-induced, temperature non- 
monotonic stick-slip friction is confirmed, prominent here 
over the (temperature- monotonic) background friction. 
The friction coefficient broadly peaks in the neighbor- 
hood of the substrate phase transition, where friction 
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FIG. 2: Friction coefficient /i = (Fx) / Fq as a function of tem- 
perature for different substrate order parameter directions A 
(correspondence as in upper panel). Hexagon occupancies in 
the lower panel illustrate the probability of the substrate "tar- 
get atom" to be found in the six valleys of the on-site potential 
after ten stick-slip events at various temperatures above and 
below Tc and for different distortion directions A. Upper and 
lower panels refer to two distinct load values. Dashed curve: 
bulk order parameter static susceptibility 



rises substantially higher than at low temperature, to 
descend again at higher temperatures beyond Tc. Sec- 
ond, there is a clear dependence of friction upon the or- 
der parameter valley in the substrate. That dependence, 
well visible for light loads at T < Tc, grows further with 
load (compare upper and lower panels in Figj2|. At large 
load and low temperature, the friction coefficient differs 
between valleys A = 1 (or 4), 2 (or 3), and A = 6 (or 5) 
by almost an order of magnitude. We stress that the fric- 
tion nonmonotony versus T is in this model of a totally 
different nature from that recently demonstrated for mul- 
tiple stick-slip regimes of motion [M] , or for multi-contact 
sliders [15]. 

So far the simulation results. Is there a linear response 
theory that may explain them? Strictly speaking the 
answer is negative, because stick-slip is a violent, non- 
linear, non-uniform perturbation. In spite of that, it is 
instructive to compare the friction simulation data with 
a different kind of linear response theory. If the slider's 
speed, jerky because of stick-slip, could be crudely re- 



placed with a large uniform speed, itself unperturbed by 
the frictional processes, then a standard "golden rule" 
linear response could be invoked, predicting an average 
dissipated power 

(P)oc^c^k|T4|V(k,Wk), (1) 

k 

where is the coefficient of the Fourier expansion of the 
slider-substrate potential, ujk = k-v, and x''(k, cjk) is the 
imaginary part of the (semi-infinite) substrate density- 
density correlation function, as in electron energy loss. 
[16]. Since k • v is a low frequency, a surge of friction 
near Tc is expected, in connection with the increased 
density of low frequency modes associated with soften- 
ing of the displacive mode, and eventually with "cen- 
tral peak" diffusive excitations in the critical regime 0. 
Without attempting to extract x''(k,co'k) from simula- 
tions, we note that since by Kramers-Kronig relations 
X'(k = 0,cj = 0) = (2/7r) ^^^^^duj, the T-dependent 
peak of at Tc, while surely not identical to that of 
x''(k, cjk), should similarly accompany the peaking dis- 
sipation in this approximation. Results in Fig. [2] in- 
deed qualitatively confirm a close similarity in the T- 
dependence of simulated stick-slip friction with that in- 
dependently obtained for the bulk susceptibility. 

Further insight in the order parameter valley depen- 
dence of friction can be obtained by inspecting the sys- 
tem's dynamics. The slider imparts mechanical kicks to 
nearby substrate atoms (details in Supplemental Mate- 
rial) following which, energy is transmitted to the sub- 
strate and degraded as Joule heat (movies provided in 
Supplemental Material). When T is low and the load 
is light the kicked substrate atoms vibrate moderately, 
harmonically, and mainly radially along the same poten- 
tial valley, see diagram (i) of Fig. |2j resulting in very low 
friction. As the load increases, still at low temperatures, 
the sliding tip causes local order parameter flips - jumps 
between valleys A ^ A' - of near-tip atoms in the sub- 
strate. The work spent by the tip to actuate this local 
flip is never returned to the tip; thus an increased flip 
rate reflects in an increased friction coefficient which is 
seen well below Tc. 

As temperature is raised, spontaneous thermal flips of 
order parameter proliferate in the substrate, eventually 
exploding critically near Tc. This is in correspondence 
with a surge of susceptibility, and to a drop of the free en- 
ergy barrier for the tip to cause additional sliding-induced 
flips; so while their number also proliferates, see Supple- 
mental Material, the friction rises to a critical maximum. 
(It should be noted however that some additional role in 
friction near Tc will be played by the muted propaga- 
tion conditions of Joule phonons injected into the sub- 
strate, where propagation may be impeded by critical 
fluctuations). Well above Tc finally all substrate atoms 
spontaneously and frequently jump over the six valleys 
(diagram (iii) Fig. |2|, offering a diminishing probabilty 
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for the slider's kicks to do work, and friction gradually 
declines. The efficiency of stick-slip in causing an or- 
der parameter flip well below Tc is clearly not the same 
for different valleys A. Fig. 5 in Supplemental Material 
shows that flips between valleys are more abundant when 
the substrate is initially polarized in A = 2 (or 1, 3, 4) 
than those with A = 6 (or 5). The force exerted by the 
slider at slip is mostly downward oriented, thus valleys 
A = 2, 3 can be kicked to A' = 1,4 (or even A' = 5,6); 
valleys A = 1,4 can be kicked to A' = 5,6; but valleys 
A = 5,6 cannot be kicked anywheres. The frictional dif- 
ferences between order parameter directions become, we 
find, even larger when the load is raised as shown in 
Fig. |2] lower panel. All observations remain essentially 
the same for a lighter mass tip, see Supplemental Ma- 
terial. Although schematic, the valley-specific efficiency 
difference in dissipation just demonstrated anticipates a 
general mechanism for the stick-slip friction dependence 
on the detailed domain orientation of the substrate or- 
der parameter, and a source of AFM frictional contrast 
between different domains. 

This observation suggests a possible use as a means to 
control friction, exemplified by the simulation of Fig. |3] 
Start out with the substrate polarized in valley A = 6, 
where the low temperature stick-slip friction is small. At 
time t = to, an external field E = [Ex{t)^ Ezit)]^ coupling 
to the order parameter in the form u • E, is turned on 
until t = ti, when it is turned off. (For a ferroelectric 
substrate, E is an electric field; for other ferrodistortive 
systems, it could be for example a uniaxial deformation). 
For sufficiently large field, the substrate overall distortion 
switches from valley A = 6 to A' = 2, and the friction 
correspondingly jumps upwards. Upon subsequent ap- 
plication of a restoring field E = [Ex(t)^ —Ez(t)\^ friction 
reverts back to low. 

In summary, we have explored the behaviour of the 
order-parameter related friction qualitatively expected 
for stick-slip sliding over a structurally 'live" substrate. 
While quantitatively model-dependent, we obtained an- 
swers to our three basic questions (i)-(iii) which by 
all signs are of wider validity. The relative magni- 
tude and detectability of the order parameter-related 
frictional effects are by necessity system-dependent and 
hard to predict. On the other hand the realm of solids 
exhibiting (nearly) continuous structural transitions is 
huge [To]. The domain contrast seen on ferroelectric Ba- 
TiO [8] and on ferroelastics such as gadolinium molyb- 
date (GMO) [H] and RbAlF4 [E] could be pursued by 
temperature studies close to the ferro-para transitions. 
Antiferrodistorters like SrTiOs (T^ - 105° K) [19^, and 
antiferroelectrics such as KMnFs (T^ - 187°i^) [20] 
would also be of interest. The choice between these or 
other materials will largely be dictated by experimental 
considerations; and so will their potential use towards 
realistic control over dry friction. 
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FIG. 3: Control of atomic stick-slip upon application of an 
external field, switching the overall distortion direction from 
A = 6 to = 2. Upper panel, applied external field versus 
time; large panel, stick-slip friction force; red and green lines, 
average friction force. Note the "braking" effect. Lower pan- 
els, magnitude and orientation of the substrate order param- 
eter. Simulations performed at ksT — 0.025[/ with Fq = 1.5 
and = 5 X 10~^. 
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